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ABSTRACT 

We have studied the nuclear activity in a sample of six quiescent early-type galaxies, with new 
Chandra data and archival fiST optical images. Their nuclear sources have X-ray luminosities ~ lO"^*- 



10^9 erg s-i {Lx/L 



Edd ^ 10~*-10~^), and colors or spectra consistent with accreting supermassive 
black holes (SMBHs) — except for the nucleus of NGC 4486B, which is softer than typical AGN spectra. 
In a few cases, the X-ray morphology of the nuclear sources shows hints of marginally extended 
structures, in addition to the surrounding diffuse thermal emission from hot gas, which is detectable 
on scales > 1 kpc. In one case (NGC 5845), a dusty disk may partially obstruct our direct view of 
the SMBH. We have estimated the temperature and density of the hot interstellar medium, which 
is one major source of fuel for the accreting SMBH; typical central densities are rig « (0.02 ± 0.01) 
cm~^. Assuming that the hot gas is captured by the SMBH at the Bondi rate, we show that the 
observed X-ray luminosities are too faint to be consistent with standard disk accretion, but brighter 
than predicted by radiatively-inefhcient solutions (e.g., advection-dominated acretion flows, ADAF). In 
total, there are ~ 20 galaxies for which SMBH mass, hot gas density, and nuclear X-ray luminosity are 
simultaneously known. In some cases, the nuclear sources are brighter than predicted by the ADAF 
model; in other cases, they are consistent or fainter. We discuss the apparent lack of correlations 
between Bondi rate and X-ray luminosity, and suggest that, in order to understand the observed 
distribution, we need to know two additional parameters: the amount of gas supplied by the stellar 
population inside the accretion radius, and the fraction (possibly ^ 1) of the total gas available that 
is accreted by the SMBH. We leave a detailed study of these issues to a subsequent paper. 

Subject headings: accretion, accretion disks — galaxies: nuclei — galaxies: individual (NGC 821, 
NGC 3377 NGC4486B, NGC 4564, NGC 4697, NGC 5845) — galaxies: structure 
— X-ray: galaxies 



1. INTRODUCTION 

There is strong evidence that all massive galaxies (not 
including dwarfs) with a spheroidal component (ellipti- 
cals, Icnticulars, and early-type spirals with bulges) host 
supermassive black holes (SMBHs) in their nuclei (Rich- 
stone et al. 1998). For a few dozen galaxies, there ex- 
ist direct determinations of their nuclear masses, either 
from measurements of spatially-resolved stellar or gas 
kinematics (for quiescent nuclei) or from "reverberation 
mapping" of their broad-line regions (for active nuclei) 
(Gebhardt et al. 2003; Ferrarese et al. 2001; Ho 2002; 
Onken et al. 2004; Peterson et al. 2004). The masses 
of SMBHs correlate well with fundamental properties of 
their host galaxies: bulge mass (Magorrian et al. 1998), 
central velocity dispersion (Ferrarese & Merritt 2000, 
Gebhardt et al. 2000), central light concentration (Gra- 
ham et al. 2001) and mass of the dark-matter halo (Fer- 
rarese 2002). However, this is not the case for the ac- 
creting luminosities of SMBHs: while a few are active 
(AGN), sometimes with bolometric luminosities close to 
or above their Eddington limits (Collin et al. 2002; Wang 
& Netzer 2003), the majority are not. X-ray observa- 
tions of these "silent" SMBHs have revealed faint sources 
with X-ray luminosities as low as ~ 10~*LEdd (Fabbiano 
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& Juda 1997; Ho 1999; Ho et al. 2001; Pellegrini et 
al. 2003a; Pellegrini 2005; Fabbiano et al. 2003, 2004). 

The bolometric luminosity can be expressed as Lboi ~ 
•qMc^, where M is the mass accretion rate and ry is the 
efficiency of conversion of gravitational energy into radia- 
tion. It is still not clear whether the extremely low lumi- 
nosity of most SMBHs is primarily due to a low accretion 
rate, or to a low radiative efficiency (i.e., 77 <C 0.1); see 
Narayan (2002) for a review. If it is due to a low accretion 
rate, we need to determine whether this is caused by a 
lack of gas available in the inner few pc, or is instead the 
result of other phenomena (for example accretion disk 
winds or jets) that may prevent the gas from accreting 
onto the SMBH. 

The most likely explanation appears to be that both 
low rj and low M play a part. It has become clear that a 
lack of fuel cannot be a general explanation: many ellip- 
tical galaxies do have large supplies of cold/ warm inter- 
stellar gas at their centers (Caon, Macchetto & Pastoriza 
2000; Macchetto et al. 1996). Fabian & Canizares (1988) 
pointed out that the accretion rate of the hot interstellar 
medium (ISM) expected from the standard Bondi (1952) 
model would still be high enough to produce X-ray lu- 
minosities at least two orders of magnitude higher than 
observed in the nuclei of most ellipticals, if accretion oc- 
curs efficiently (77 ^ 0.1). Recent high-resolution Chan- 
dra observations have led to the measurement of both the 
nuclear X-ray emission and of the physical parameters of 
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the surrounding hot ISM, greatly strengthening this con- 
clusion (e.g., Pellegrini 2005 and references therein). 

Advection of energy into the SMBH (e.g., Narayan & 
Yi 1994, 1995) is a plausible mechanism for low radia- 
tive efficiency, and may be the explanation for many ex- 
tremely sub-Eddington sources. However, it cannot be 
the only explanation for all sources, either. For example, 
advective accretion cannot be invoked to explain some 
very under-luminous stellar-mass compact objects in our 
Galaxy (Perna ct al. 2003). Low Al in addition to low 
?7 has also been suggested as the major reason for the 
overall faintness of Sgr A* (Bower et al. 2003; Quataert 
2003). 

Hence, there is now considerable observational evi- 
dence that quiescent or very faint SMBHs are at least 
partly due to low radiative efficiency, and partly to accre- 
tion rates onto the compact object lower than predicted 
from the amount of gas available. However, the physical 
processes responsible for this reduced accretion rate are 
still not well understood, and various competing models 
have been proposed. Hydrodynamical simulations (e.g., 
Perna et al. 2003, their Sect. 2.3 and references therein) 
have shown that disk winds and outflows may prevent 
most of the inflowing gas from reaching the compact ob- 
ject, so that only a small fraction of the available mass is 
accreted. Another possibility is that convection sets up a 
barrier to accretion, by increasing the pressure and pre- 
venting the gas from falling in supersonically (Narayan, 
Igumenshchev, & Abramowicz 2000). Alternatively, Tan 
& Blackman (2005) speculated that efficient cooling in- 
side the Bondi radius causes the infalling gas to condense 
into stars in a gravitationally unstable disk, instead of 
sinking into the SMBH. 

Another promising scenario is that most of the accre- 
tion power is not radiated, but is instead carried away 
by a steady, relativistic jet, with bulk kinetic luminos- 
ity Lj - O.lAfc^ > rjMc^. In some cases (e.g., M87: 
di Matteo et al. 2003, IC 4296: Pellegrini et al. 2003b; 
IC 1459: Fabbiano et al. 2003), it has indeed been shown 
that including the kinetic energy of the jet reconciles 
the observed energy budget with the expected accretion 
rate. If a relativistic jet is present, synchrotron radio 
emission is also expected. The "fundamental-plane" re- 
lation (Merloni, Heinz, & Di Matteo 2003) between the 
observed radio luminosities. X-ray luminosities, and BH 
masses supports the hypothesis that inactive SMBHs 
have low accretion rates as well as low radiative effi- 
ciency (?7 < 10^"^), and that most of the accretion power 
is carried outwards by a steady jet. However, in other 
cases (e.g., NGC821: Fabbiano et al. 2004) it is more 
likely that the accretion power is released in intermit- 
tent, short-lived outbursts (Binney & Tabor 1995; Ciotti 
& Ostriker 2001; Jones et al. 2002; Omma et al. 2004); if 
that is the case, the relative fraction of SMBHs observed 
in an active or quiescent state would depend on the du- 
ration of the duty cycle. Alternatively, in a few cases 
the X-ray emission from the SMBH may be absorbed by 
dust (often found in the inner ^ 100 pc of elliptical galax- 
ies: Rest et al. 2001) and re-emitted in the infrared (IR). 
Finally, it is still possible that a few luminous ellipti- 
cal galaxies with partially depleted stellar cores may not 
have large gas supplies to feed the SMBH. Determining 
which of these physical processes are at play in quiescent 



SMBHs will also help us understand and quantify the 
SMBH feedback onto its host galaxy. 

For a sample of 18 nearby galaxies with nuclear X-ray 
luminosity Lx and Bondi mass accretion rate of hot gas 
estimated from Chandra data, Pellegrini (2005) found 
that Lx can vary by orders of magnitude and shows a 
large scatter when plotted against BH mass or Bondi 
rate. While in some cases the observed luminosities are 
consistent with the prediction of advective models, for 
other nuclear sources the observed luminosities are much 
lower, and therefore scenarios including outflows or con- 
vection are needed. A scenario with feedback from ac- 
cretion on the surrounding gas was considered the most 
promising. 

1.1. Objectives and targets of our study 

To test these scenarios, we chose a further sample of 
six nearby, quiescent elliptical galaxies with direct SMBH 
mass determinations, from Gebhardt et al. (2003) and 
Kormendy et al. (1997). Stringent upper Hmits on AGN 
activity exist for all these SMBHs, from radio, optical 
(Balmcr emission lines), and previous X-ray studies (Ho 
2002), and from a lack of UV bumps. Our goal is to 
explore the most extreme cases of sub-Eddington nu- 
clear activity, so we only focused on galaxies with SMBH 
masses > lO^Af0 and X-ray luminosity < 10^° erg . 
We also excluded galaxies more distant than 30 Mpc. 

In five of the six galaxies, the SMBH mass has been 
accurately determined with three-integral models (Geb- 
hardt et al. 2003)^. The only exception is NGC4486B, 
whose mass estimate was based on two-integral models 
(Kormendy ct al. 1997) and is therefore much more mi- 
certain. For example, a revised error range in Kormendy 
& Gebhardt (2001) allows for masses between 2 x 10^ 
and IQ^Mq. A mass w 7 x lO^M© was quoted as a pre- 
liminary result of an ongoing study by Gebhardt (2005, 
priv. comm.). As a corollary of our study, we have ob- 
tained an alternative (indirect) SMBH mass determina- 
tion for NGC4486B (Soria et al. 2005 = Paper II), using 
the MBH-velocity dispersion (Tremaine et al. 2002) and 
MBH^Sersic index (Graham et al. 2003, 2005 in prep.) re- 
lations. We are aware that the larger mass uncertainty, 
and its pecuhar morphological type, make NGC4486B 
somewhat different from the rest of our Chandra targets, 
among which it was included for historical reasons. Clas- 
sified as a compact elliptical (cE) galaxy, NGC4486B 
may have originated as the bulge of a much larger, tidally 
stripped disk galaxy (Jerjen & Binggeli 1997; see also 
Graham 2002). 

Our main objective is to determine the energy and 
mass budget for the nuclear regions of our sample galax- 
ies, and for a few others for which SMBH mass and lu- 
minosity are available in the literature. We want to an- 
swer three fundamental questions: how much gas is avail- 
able for accretion? how much energy is radiated by the 
SMBH? where does the rest of the mass and energy go? 
Moreover, this may also help to explain why most of the 
SMBHs in the local universe are very faint or quiescent. 

In three of those five nuclei, the dynamical sphere of influence 
of the SMBH has been resolved by HST (Mcrritt & Fcrrarese 2001; 
Ho 2002), making the mass determination more reliable; in the re- 
maining two cases, NGC 821 and NGC 3377, the sphere of influence 
is not resolved. 
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We begin by using Chandra data to investigate the spa- 
tial and spectral properties of the nuclear X-ray emis- 
sion. By comparing the observed X-ray luminosity of 
the SMBH with the temperature and density of the sur- 
rounding gas, we can constrain some of the accretion 
parameters. In particular, if we could precisely deter- 
mine the mass accretion rate onto the SMBH, we would 
find the efficiency. However, a direct measurement of 
the accretion rate is still beyond today's observational 
capabilities, given the characteristic size of SMBH event 
horizons (~ 10~^ pc). So, we need to turn the problem 
around. We shall use model predictions for the radiative 
and non-radiative efficiency (e.g., standard disks, or ad- 
vective accretion) and compare them with the observed 
luminosity and gas density available inside the SMBH 
accretion radii (characteristic size ^ 10 pc). From this, 
we will infer what fraction of the gas available is really 
accreted. In Paper I, we shall use Chandra to estimate 
the hot gas component to the accretion flow; we shall 
look for correlations between the classical Bondi inflow 
rate from the hot ISM, and the observed SMBH luminos- 
ity. In Paper II, we shall use optical data to determine 
whether additional sources of fuel (e.g., warm gas from 
stellar winds) may have been missed by the X-ray study. 
We shall then discuss the correlation between the total 
gas accretion rate and the observed SMBH luminosities. 
Finally, we shall discuss how to balance the mass and 
energy budget, taking into account total gas injection, 
accretion onto the SMBH, and outflows. 

2. X-RAY OBSERVATIONS AND DATA ANALYSIS 

The six elliptical galaxies that are the object of our 
study (Table 1) were observed with Chandra (Weisskopf 
et al. 2000) ACIS-S, centered on chip 7 (i.e., ACIS- 
S3). For one of the galaxies, NGC4697, we used all 4 
observations currently available in the public archive. 
We analysed the data with the CIAO v3.2 software^ 
for NGC4697, and with CIAO v3.1 for the other tar- 
gets; we checked that the difference between the two 
CIAO versions is negligible for the purposes of this 
work. We downloaded the bias files for each observa- 
tions, and built new bad-pixel files^. We then built 
new level=l and level=2 event files with the CIAO script 
acis_process_everLts, which includes the application of 
a time-dependent gain correction to the data. (The 
time-dependent gain adjustment was not applied to the 
2000 observation of NGC 4697, because it was taken 
with a focal plane temperature different from the nor- 
mal — 120 C). We also improved the astrometric accu- 
racy by calculating and removing the (negligible) aspect 
offsets, using the CIAO Aspect Calculator''. No signifi- 
cant background flares were observed in these data, so 
no time filtering was necessary. The background level 
was generally very low except for the observation of 
NGC4486B. The energy range was filtered to 0.3-10 
keV, with the CIAO task dmcopy. We identified the point 
sources with the wavdetect task, and for the brightest of 
them, we extracted their spectra with psextract (which 
also builds response and auxiliary response files, with 
the mkrmf /mkarf tasks). After subtracting the discrete 

^ http://cxc.harvard.edu/ciao/ 

® Thread at http:/ / cxc.harvard^edu/ciao/threads/acishotpixels^__^ 
^ Thread at http: / / cxc.harvard^Su/cmo/threads/arcseccorrectionTI 



sources, we used the CIAO script acisspec to extract 
the spectra of diffuse emission from extended regions, 
together with their weighted response files and back- 
grounds. We then used XSPEC 11.3.1 (Arnaud 1996) for 
the spectral fitting of point- like and extended sources*. 

3. RESULTS 

We focus in particular on the new results for 
NGC 3377, NGC4486B, NGC 4564 and NGC 5845. The 
remaining two galaxies have already been studied in de- 
tail (NGC 821: Fabbiano et al. 2004; NGC 4697: Sarazin 
et al. 2001): here we re- analyse and discuss their nu- 
clear X-ray properties in the context of our new sample 
of massive, quiescent elliptical galaxies. In particular, 
for NGC 4697 we use both the 2000 Chandra observation 
that was studied by Sarazin et al. (2001), and three more 
recent observations that became public in early 2005. 

3.1. Morphology of the nuclear X-ray sources 

In all six galaxies we detect X-ray emission at and 
around the nuclear position. In five cases (all except 
NGC 821, see Fabbiano et al. 2004), the centroid of 
the nuclear X-ray source coincides with the optical/IR 
nucleus identified in the Two-Micron All Sky Survey 
(2MASS) and United States Naval Observatory (USNO) 
Catalogues, within the astrometric accuracy of Chandra 
(the 90% uncertainty circle has a radius of < 0'.'6). For 
NGC 4697, relative astrometry with an accuracy better 
than w 0'.'3 can be obtained, because a few non-nuclear 
X-ray sources have point-like optical counterparts (glob- 
ular clusters). 

Firstly, we analyzed the morphology of the nuclear 
X-ray sources to determine whether they are consistent 
with point-like emission from the SMBH or are instead 
extended. From the theoretical ACIS-S3 point spread 
function and from our fitting of other point sources de- 
tected in the central region of the chip, we expect a Gaus- 
sian FWHM = 0'.'89 ± 0'.'15 for a point-like source in the 
0.3-8 keV energy band. The (weak) nuclear source of 
NGC 4486B is consistent with being point-hke (FWHM 
w 0'.'9), implying a size < 70 pc; it is also well isolated: 
no other X-ray sources are detected in that galaxy. The 
nuclear source in NGC 4697 (Figure 1, bottom left) is 
also point-like (size < 50 pc) but is located in a region 
with high source density, with at least another 5 point- 
like sources detected within a radius of 300 pc. Their 
spatial distribution around the nuclear source is proba- 
bly a chance grouping of X-ray binaries in the galactic 
bulge. 

The nuclear source of NGC 4564 (Figure 1, top right) 
has a FWHM « 2'.'2 « 160 pc in the 0.3-8 keV band; 
therefore, it is probably extended. In the remaining three 
galaxies, NGC 3377, NGC 5845 and NGC 821, struc- 
ture in the nuclear region is more clearly visible. In 
NGC 5845 (Figure 1, bottom right), the nuclear emis- 
sion has an irregular, "fuzzy" appearance over a region 
« 600 pc across. In the other two galaxies, NGC 821 and 
NGC 3377, the emission has instead an elongated, "jet- 
hke" appearance. In NGC 3377 (Figure 1, top left), an 
elongated feature seems to link the nuclear source with 
another, fainter source located « 3" « 160 pc to the 
East. For a discussion of the elongated substructure in 

* http:/ /heasarc. gsfc.nasa.gov/lheasoft/xanadu/xspec/ 
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TABLE 1 
Targets of our study 



Galaxy 
(1) 


Type 
(2) 


d 

(Mpc) 
(3) 


Mb 
(mag) 
(4) 


(1O»M0) 
(5) 


Obs. Date 
(6) 


Obs. ID 
(7) 


Exp. time 

(ks) 
(8) 


NGC 821 


E6 


24.1 ±2.0 


-20.7 


gc:+0.36 
>J-O3_0.35 


2002/11/26 
2002/12/01 


4408 
4006 


25.3 
13.7 


NGC 3377 


E5-6 


11.2 ±0.5 


-19.2 


i.oi?,? 


2003/01/06 


2934 


39.6 


NGC 4486B 


cEO 


16.9 ± 1.3 


-16.8 


[6.oll:°] 

5+°-^ 
U.O_Q 2 


2003/11/21 


4007 


36.2 


NGC 4564 


E6/S0 


15.0 ± 1.2 


-19.0 


U.OO_g Qg 


2003/11/21 


4008 


17.9 


NGC 4697 


E6 


11.7 ±0.8 


-20.3 


1 7+0.2 
^•'-0.1 


2000/01/15 
2003/12/26 
2004/01/06 
2004/02/12 


784 
4727 
4728 
4729 


39.3 
39.9 
35.7 
38.1 


NGC 5845 


E3 


25.9 ±2.7 


-18.8 


2 4+° '* 


2003/01/03 


4009 


30.0 



Note. — Col.(l): see Fabbiano et al. (2004) for a detailed study of NGC 821, and Sarazin, Irwin & Bregman (2001) for the analysis of the first 
observation of NGC4697. Col. (2): from the NASA/IPAC Extragalactic Database (NED). Col. (3): from Tonry et al. (2001), except for NGC4486B, 
from Neilsen & Tsvetanov (2000). Col.(4): from NED, adopting the distances in Col. (3). Col. (5): from Gebhardt et al. (2003), except for NGC 821 
(Richstone et al. 2004), and for NGC4486B. The latter has two alternative values. The higher mass (in square brackets) is from Kormendy et 
al. (1997), using a two-integral model. A similar result (5.0^^'g X 10^ Mq) was obtained with the same method by Kormendy & Gebhardt (2001). 
However, the reliability of this value is questionable, both because of the method on which it is based, and because it is uniquely discrepant from 
the BH mass inferred from all indirect indicators (i.e., the global correlations with stellar velocity dispersion, Scrsic index, bulge mass). The 
alternative, indirect mass determination was obtained in our Paper II (Soria et al. 2005), from the global correlation between Sersic index in the 
optical brightness profile and BH mass (Graham et al. 2003, 2005 in prep.). 



TABLE 2 

ACIS-S NET COUNTS AND COUNT RATES FROM THE NUCLEAR SOURCES 



Galaxy 


0.3-10 keV 


ACIS-S Net Counts 
0.3-1 keV 1-2 keV 


2-10 keV 


ACIS-S Count Rates (10""' s" 
0.3-10 keV 0.3-1 keV 1-2 keV 


2-10 keV 


NGC 821 (core) 








not detected 
















NGC 821 ("jet") 


137.2 ± 26.4 


34.3 ±7.7 


63.5 ±9.5 


39.4 ± 23.4 


36.2 ± 7.0 


9.0 


±2.0 


16.8 


±2.5 


10.4 


± 


6.2 


NGC 3377 


113.1 ±11.6 


41.5 ±6.8 


45.5 ±6.8 


26.1 ±5.4 


29.3 ±3.1 


10.9 


±1.8 


11.5 


±2.4 


6.9 


± 


1.4 


NGC4486B 


38.2 ±8.2 


27.7 ± 6.2 


6.8 ± 7.0 


3.7±5.1 


10.6 ±2.3 


7.7 


±1.7 


1.9 


±1.9 


1.0 


± 


1.4 


NGC 4564 


72.3 ±9.3 


28.3 ±5.9 


27.5 ±5.6 


16.5 ±4.4 


43.7 ±5.6 


17.0 


±3.5 


16.7 


±3.4 


10.0 


± 


2.7 


NGC 4697 (1) 


128.9 ± 11.6 


57.3 ±7.8 


40.4 ±6.4 


31.4 ±5.7 


32.8 ±3.0 


14.6 


±2.0 


10.3 


± 1.6 


8.0 


± 


1.4 


NGC 4697 (2) 


288.0 ± 17.3 


97.0 ± 10.2 


120.7 ± 11.1 


69.3 ±8.5 


25.3 ± 1.5 


8.5 


±0.9 


10.6 


± 1.0 


6.1 


± 


0.7 


NGC 5845 


117.8 ± 11.8 


36.5 ± 6.8 


56.9 ± 7.9 


24.4 ±5.3 


39.4 ±3.9 


12.2 


±2.3 


19.0 


±2.7 


8.2 


± 


1.8 



XoTE. — The values for NGC 821 are from the combined 2002 Nov/Dcc observations; sec Fabbiano et al. (2004) for details on the definition of 
the jet-like nuclear feature as opposed to a point-like (non-detected) core. The first line for NGC 4697 is from the 2000 Jan observation, and the 
second is for the combined 2003 Dec-2004 Feb observations; the reduction in count rate is due to the degraded performance of ACIS-S. 



the nucleus of NGC 821, see Fabbiano et al. (2004). Lim- 
iting our imaging study to the 1.5-8 keV band (Figure 
2), we still find some extended structure in NGC 5845. 

3.2. X-ray colors and spectra of the nuclear sources 

We extracted X-ray spectra of the nuclear sources, us- 
ing the extraction regions plotted in Figure 1 (dashed 
white circles); the background regions were chosen as an- 
nul! around the nuclear sources, and wc took care of ex- 
cluding nearby sources. For NGC 4697, we extracted sep- 
arate spectra for the four observations. For some tasks, 
it is possible and convenient to coadd the 2003 Dec-2004 
Feb observations and study them together; however, the 
analysis of the first dataset (2000 Jan) has to be carried 



out separately because of the ACIS-S sensitivity degrada- 
tion, particularly at lower energies, over those four years. 

Firstly, we compared the broad-band X-ray colors (Ta- 
ble 2), before introducing any spectral modelling. In five 
of the six nuclear sources, the X-ray colors are consistent 
with typical AGN spectra, a power-law spectrum of index 
« 1.5-2 (Figure 3). The only exception is NGC4486B, 
which has softer colors, with most of the counts detected 
below 1 kcV, although the errors are large because of the 
relatively high background level and the small number of 
counts. 

We then estimated the emitted fluxes and luminosities 
by assuming an absorption column density equal to the 
Galactic line-of-sight (Table 3; from Dickey & Lockman 
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1990). The color-color diagram (Figure 3) and the X- 
ray spectra (at least for the nuclear sources with enough 
counts) justify this assumption, suggesting that the to- 
tal column density is not much higher than the Galactic 
value; hence, we are only slightly underestimating the 
emitted luminosities. We fitted the (ungrouped) spectra 
with an absorbed power-law (Table 3), using the Cash 
statistic (Cash 1979)^. We also grouped the same spec- 
tra to achieve a minimum signal-to-noise ratio of 3 in 
each spectral bin, and fitted them with the statistics 
(Table 3). For sources with a very limited number of 
counts, we consider the results based on the Cash statis- 
tic to be more reliable. For NGC 4697, we increased the 
signal-to-noise ratio by coadding the spectra from 2000 
and 2003-2004, using the algorithm of Page, Davis & 
Salvi (2003), which generates a count- weighted average 
of their response matrices (Figure 4). From a power- 
law fit (Table 3) we calculate an emitted luminosity 
Lx = (3.9 ± 0.2) X 1038 erg s'^ in the 0.3-10 keV band; 
we also checked that there are no significant spectral or 
luminosity changes between the two epochs. The emit- 
ted luminosities of the nuclear sources in the other five 
galaxies are between ?» 3 x 10^^ and « 3 x 10"^^ erg s~^, 
in the same band (Table 3). 

For NGC4486B, we have also tried fitting the nu- 
clear source with a thermal-plasma model, to account 
for the possibility that it is a very compact patch of 
thermal emission from the ISM. We obtain a best-fit 
kT = 3.7tl°Q keV, but the fit is statistically worse (C- 
statistics = 54/297) than for the power-law model (C- 
statistics = 46/297). Henceforth, we shall assume for 
this work that most of the X-ray emission comes from 
the SMBH. 

3.3. Scattering and obscuration of the nuclear SMBH? 

The nuclear X-ray source in NGC 5845 appears ex- 
tended on the scale of a few hundred pc, larger than the 
FWHM of the Chandra/ ACIS PSF for a point source 
(Section 3.1). The count rate is insufficient to map the 
X-ray colors or the spectral properties across the nuclear 
region. However, we can at least say that the "fuzzy" 
emission around the nuclear position is not significantly 
softer than a typical AGN-like power-law spectrum: the 
emission is extended even in the hard band (2-10 keV). 
The simplest explanation is that the emission is due to 
a few faint, unresolved LMXBs in the inner region. If 
that was the case, the true luminosity of the SMBH in 
NGC 5845 would be a few times lower than estimated in 
Table 3. 

Alternatively, we note that the X-ray morphology is 
reminiscent of the extended X-ray emission associated 
with obscured Seyfert nuclei (for NGC 4151: Elvis, Briel 
& Henry 1983; Elvis et al. 1990; Ogle et al. 2000; for 
NGC 1068: Ogle et al. 2003), which is understood as 
scattered nuclear emission. Dust is known to be present 
in the nuclear region of most ellipticals, especially those 
with disky isophotes (van den Bosch et al. 1994; Kor- 
mendy et al. 2005 and references therein). The nuclear 
region of NGC 5845 contains a dusty disk of radius w 100 

^ This statistic can be used in the same way as the statistic 
to find confidence intervals. One finds the parameter values that 
give C = Cjnin + N, where N is the same number that gives the 
required confidence for the number of interesting parameters as for 
the case (Arnaud 1996). See the XSPEC manual for details. 



pc (Quillen, Bower & Stritzinger 2000), i.e. an order of 
magnitude larger than the Bondi accretion radius. In 
fact, NGC 5845 is one of only two galaxies known for 
containing both a dust disk and an associated, nuclear 
stellar disk (the other one is NGC 4486 A). It has been 
suggested (Kormendy et al. 2005; Kormendy et al. 1994) 
that, as cold gas and dust settle towards the center, they 
form stars and build a stellar disk. 

In NGC 5845, the disk is seen almost edge-on (Figure 
5). Thus, it is possible that most of the direct X-ray emis- 
sion from the SMBH along our line of sight is obscured 
by dust, and most of the observed, extended emission 
comes from the surrounding photo-ionized plasma which 
is scattering the nuclear continuum. If this was the case, 
the true SMBH luminosity in NGC 5845 could instead 
be at least an order of magnitude higher than estimated 
in Table 3. Given the apparent thin-disk geometry of 
the dusty feature (Figure 5), we may expect that the 
disk/torus can intercept and re-radiate in the far-infrared 
a fraction ~ 10-20% of the bolometric luminosity. If the 
far- infrared nuclear luminosity > 10**^ erg s^^, it is po- 
tentially detectable by Spitzer (flux density ^ 10 mjy in 
the Spitzer/MlPS bands). This will be investigated in a 
follow-up work, for which Spitzer time has already been 
obtained. For the present paper, we assume that the ob- 
served X-ray luminosity of the SMBH in NGC 5845 is a 
good approximation of the emitted luminosity. 

A dusty disk or torus of outer radius 200 pc is also 
apparent in the nuclear region of NGC 4697 (Figure 6). 
However, obscuration or scattering of the central X-ray 
source may not be as significant as in NGC 5845 because 
our viewing angle is ~ 75°. A more detailed investigation 
of the physical processes at work in a dusty nuclear disk 
is beyond the scope of this paper. 

3.4. Jet-like features 

The morphologies of the X-ray emission in the nuclear 
regions suggest possible jet-like features in two of our tar- 
get galaxies, NGC 821 and NGC 3377. The former was 
discussed in Fabbiano et al. (2004), who suggest that 
the X-ray flux is consistent with synchrotron emission in 
a jet, or with hot thermal plasma shocked by intermit- 
tent nuclear outbursts. For NGC 3377, we compared the 
Chandra image with an archival _ffS'TyWFPC2 image in 
the F555W filter, to determine whether the elongated 
feature in the X-ray image has an optical counterpart; 
however, we do not see any obvious distortions in the 
optical isophotes (Figure 7, top panel). 

We then searched for possible fainter optical signatures 
in the nuclear region of NGC 3377. We determined the 
mean 2-D surface brightness distribution of the central 
region of NGC 3377 with the help of a boxcar pixel- 
smoothing technique and an isophotal fitting routine 
written in IRAF. Both techniques were recently employed 
to uncover spiral and bar structures in seemingly nor- 
mal dwarf elliptical galaxies in the Virgo (Jerjen, Kalnajs 
& Binggeli 2000; Barazza, Binggeli & .Jerjen 2002) and 
Coma clusters (Graham, Jerjen & Guzman 2003). The 
reeonstrticted model galaxy is subtracted from the orig- 
inal, leaving just the small-scale substructures in the 
residual image. The result for NGC 3377 derived in this 
manner is shown in the bottom panel of Figure 7: we 
find a small, brighter stellar disk but no feature corre- 
sponding to the X-ray elongation, neither parallel not 
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Fig. 1. — ACIS-S image (0.3—8 keV band) of the nuclear regions in four of the galaxies showing extended structures. (For the other two 
galaxies in our sample: see Fabbiano et al. (2004) for NGC821; the nuclear source is point-like in NGC4486B). In each panel, the dashed 
white circle is the source region used to determine the nuclear luminosity, spectrum and colors (Tables 2, 3 and Figure 3). The green circles 
(radius of 0'.'6) are the 2MASS positions of the galactic nuclei, except for NGC4697, where, for consistency with previous work, the green 
circle marks the average between the 2MASS and the USNO Bl.O Catalogue positions (the difference between the two is 0'.'4). North is 
up. East to the left. 



perpendicular to the disk plane. 

None of the six galaxies in our sample have signifi- 
cant radio detections; however, five have at least reliable 
upper limits to their radio-core emission, from VLA ob- 
servations. In NGC821, the radio-core flux at 5 GHz 
is < 0.3 mJy (Wrobel & Heeschen 1991); in NGC3377, 
NGC4697 and NGC5845, the radio flux at 8.4 GHz is 
< 0.1 mJy (Krajnovic & Jaffe 2002). In NGC4564, the 
upper limit is 0.52 mJy at 8.4 GHz (Wrobel & Herrn- 
stein 2000); however, the core is marginally detected at 
the level of « 2 mJy at 1.4 GHz, in the NRAO/VLA Sky 
Survey. A comparison of radio and X-ray nuclear fluxes 
is discussed in Paper II. 

4. X-RAY EMITTING, DIFFUSE GAS 

4.1. X-ray brightness profiles 

Unresolved X-ray emission in elliptical galaxies is a 
combination of faint low-mass X-ray binaries (LMXBs) 
and hot gas, with kT - 0.3-0.7 keV (Trinchieri & Fab- 
biano 1985; Blanton, Sarazin & Irwin 2001; Sivakoff, 
Sarazin & Irwin 2003; Randall, Sarazin & Irwin 2004). 
In this work, we focus on the nuclear BH activity and 
do not set out to study the X-ray source population and 
large-scale structure of the galaxies. However, we still 



need to investigate the diffuse hot gas component, be- 
cause it provides a constraint to the energy budget of 
the nuclear BH. From the surface brightness and temper- 
ature of the diffuse emission, we can in principle estimate 
the density of the hot gas in the nuclear region; hence, 
wc can obtain a lower limit to the gas density available 
for accretion onto the BH, which can be compared to the 
observed luminosity, to constrain the accretion efficiency. 

To obtain the surface brightness of the unresolved 
emission (and, in particular, the thermal gas compo- 
nent), firstly we filtered the event files to remove the 
point sources detected with wavdetect. In the case of 
NGC4697, we combined all four observations to reach 
a deeper detection limit for the point sources; then we 
analysed the unresolved emission separately for the 2000 
and 2003-2004 datasets. For each galaxy, we created a 
stack of concentric annular regions extending up to 120" 
from the nucleus, and a background annulus outside this 
region. We used the CIAO task dmextract to obtain the 
background-subtracted radial profiles. Since we are only 
interested here in the diffuse gas contribution, we re- 
stricted the energy range to 0.3-1.5 keV, where the hot 
thermal plasma component is relatively stronger, com- 
pared with the unresolved LMXB contributions. Above 
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Fig. 2. — As in Figure 1, but limited to the 1.5—8 keV band. The green circles (radius of 0'.'6) mark the best available optical/IR position 
for the nuclei. North is up, East to the left. 



TABLE 3 

Spectral parameters and emitted luminosities from the nuclear sources. All spectra have been fitted with an absorbed power law, and 
column density fixed to the line-of-sight galactic value. for each source, the first fit is based on the cash statistics, the second (when 

applicable) on THE STATISTICS. ERRORS ARE 90% CONFIDENCE LEVELS FOR ONE INTERESTING PARAMETER (Ax^ = 2.7) 



Galaxy 


Galactic 
(1020 cm-2) 


Photon index F 


Normalization at 1 keV 
(10"'' ph. keV-i cm-2 s~i) 


X^/dof 


.tj0.3-10kcV 
(10^8 erg s-i) 


NGC821 ("jet") 


6.4 


^■^'-'-0.31 


''■'-'-1.0 




25.5 ±2.3 






1 5n+°-^i 

1.0U_g 29 


4.9 ± 1.0 


7.5/6 


24.4 ±2.2 


NGC 3377 


2.9 




O q + 0.6 

''•''-0.6 




3 fi+0-7 
''•"-0.8 






1.8511:^ 


''•^-1.8 


1.6/7 


o 4+2.6 


NGC 44863'= 


2.6 


2 4+°-^ 


1 6+0-5 






NGC 4564 


2.4 


2 4+°-'* 


4 5+1-1 
4.0_j 2 










1 6+1-5 
i.0_2 1 


''•''-3.7 


0.76/4 


I0.5l25g3 


NGC 4697^* 


2.1 


, Ko+0.23 
1.00_Q 22 


2.9«-J 


6.6/22 


o q+0.2 
j.y_Q 2 


NGC 5845 


4.3 


2 n3+°-2* 

Z.UO_Q 29 






24.0tli 






1 4+0.8 


5 1+2-1 

^•1^-2.2 


1.3/7 





^from Dickey & Lockman 1990. 

^Emitted luminosity inferred from the spectral fit. 
■^Fitted with Cash statistics only. 
■^Fitted with the statistics. 



~ 1.5 keV, the hot gas component is, mstead, neghgible, 



and the unresolved emission is dominated by the power- 
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-0.6 -0.4 -0.2 0.2 

(H-M)/(S+M + H) 

Fig. 3. — X-ray color-color diagram for the nuclear emission 
in the six galaxies of our sample. We overplotted the expected 
colors for an optically-thin thermal plasma model and for power- 
law models of different photon indices, with varying line-of-sight 
column density (ran running from 2 X 10^0 to 5 X lO*' cm along 
each curve, from the bottom to the top). The galaxies are labelled 
from 1 to 6 as in Table 1. For NGC821, for which no point-like 
source is detected at the nuclear position, we plotted the average 
colors of the elongated, jet-like feature. For NGC4697, we only 
used the 2003-2004 data, to minimize the effect of color changes 
due to the sensitivity degradation. The model curves are the colors 
expected for an observation in mid-2003. 




0.5 1 2 5 10 

channel energy (keV) 



Fig. 4. — Spectrum of the nuclear source in NGC4697, with 
best-fit model and residuals. The data have been fitted with an 
absorbed power law, of photon index F = 1.5 it 0.2 (Table 3). 



law LMXB component. 

To determine the physical properties of the extended 
X-ray emission, we need to measure the radial profiles 
of the gas surface brightness S(r), and the central sur- 
face brightness (a distance-invariant quantity). The 
surface brightness S{t) does not come directly from the 
observations: what we have instead is the surface density 
S'(r) of the count rates in the 0.3-1.5 kev band, for the 
total unresolved emission, including contributions from 
both the interstellar gas and faint LMXBs. We show 
S"(r) for five of the six galaxies in our sample (Figures 8, 
9). The strongest emission is found in NGC4697 (Fig- 
ure 10; see also Sarazin et al. 2001). On the other end, no 
unresolved emission is detected in NGC4486B, though 
the relatively higher background level of that observa- 
tion makes it more difficult to obtain a tight upper limit. 




Fig. 5.— Greyscale i/ST/WFPC2 image (F555W filter) of the 
inner region of NGC 5845, with X-ray contours of the nuclear emis- 
sion overplotted (cf. Figure 1, bottom right panel). North is up. 
East to the left. The nucleus contains an 100-pc-radius dusty 
and stellar disk, which may contribute to obscuring or scattering 
the central X-ray source. In an alternative scenario, the dusty disk 
may reduce the accretion rate onto the BH by storing cold gas and 
triggering star formation. 



To calculate the density of the hot gas, we need to con- 
vert the observed count rates into physical flux units, and 
to estimate what fraction of the unresolved emission is 
due to diffuse gas. This requires a spectral analysis. 

Only two galaxies (NGC 4697 and NGC 3377) have 
enough detected counts in the unresolved X-ray emission 
to allow for spectral fitting. After excluding the discrete 
sources found by wavdetect, we used the CI AO script 
acisspec to extract the spectrum of the unresolved emis- 
sion from a 90" circle, taking an annulus between 120" 
and 150" as background region. For NGC 4697, we ex- 
tracted separate spectra (with their respective response 
files) for the 2000 and 2003-2004 observations; we then 
combined the two spectra using the algorithm of Page 
et al. (2003). We fitted the spectra in XSPEC, assuming 
an absorbed power-law plus a single-temperature ther- 
mal plasma model, wabs*(vapec + po). First, we fixed 
the absorption column density to the Galactic value, the 
metal abundance of all elements to the solar value, and 
the photon index of the power-law component to 1.7 (av- 
erage value for LMXBs), to reduce the number of free 
parameters (Tables 4, 5). This already produces a good 
fit (x^ = 43.0/51) As previously noted {ROSAT study: 
Pellegrini & Fabbiano 1994; Chandra study: Sarazin 
et al. 2001), the interstellar gas in NGC 4697 is cooler 
[kT K, 0.3 keV) than is usually found in X-ray bright el- 
lipticals; this may be the consequence of a shallow grav- 
itational potential. We further improved the fit by vary- 
ing the metal abundances of the main elements: we find 
that an overabundance of C and N over the a elements 
by a factor of ~ 5-10 significantly reduces the below 
1 keV. No significant Mg or Si lines are detected (Fig- 
ure 11), an indication that the X-ray emitting gas may 
not be strongly metal-enriched from Type-II SNe. The 
gas in NGC 3377 is consistent with a higher temperature 
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Fig. 6.— Greyscale HST/ ACS image (F475W filter) of the inner 
region of NGC 4697 (top panel), with a close-up view of the nuclear 
region (bottom panel). In both panels, North is up, East to the 
left. The X-ray contours are on a square-root scale with arbitrary 
zeropoint. We improved the relative astrometry of the HST and 
Chandra images by taking advantage of a few X-ray sources with 
an optical counterpart (possibly globular clusters in NGC 4697). 
Four of them are visible here (in addition to the nuclear source), 
labelled as "1" through "4". The optical nucleus is point-like in 
HST (size < O'.'l ~ 5 pc) and is probably a stellar cluster. 



{kT w 0.6 keV), although the uncertainty is larger, due 
to the limited number of detected counts (only « 130; 
Figure 12). There are not enough counts to vary the 
metal abundances. 

By fitting a spectral model to the unresolved X-ray 
emission, we can estimate the conversion factors between 
the observed count-rate density S'{r) (units of counts 
s~^ arcsec"^ in a specific band), and the bolometric 
surface brightness S{r) of the diffuse hot gas. For the 
spectral model fitted to NGC 3377, we obtain a conver- 
sion factor of 1 count arcsec^^ in the 0.3-1.5 keV 
band — (2.5 ± 0.5) erg cm~^ s~^ arcsec"^; ~ 50% of 
the unresolved counts in the 0.3-1.5 keV band are from 
truly diffuse emission. This is an average value over the 
whole extraction region; in fact, the diffuse emission may 
be more centrally peaked (Kim & Fabbiano 2003). For 
the 2000 observation of NGC 4697, 1 count s"^ arcsec"^ 
= (4.8 ± 0.4) erg cm-^ s'^ arcsec-^; in 2003-2004, 1 




Fig. 7. — Comparison between X-ray (Chandra/ ACIS, 0.3—8 keV 
band) and optical (_ff5'T/WFPC2, in the F555W filter) images of 
the inner region of NGC 3377. Top panel: the elongated X-ray 
emission feature at the nuclear position does not have any corre- 
spondence in the optical isophotes. Bottom panel: X-ray emission 
contours overplotted onto a greyscale residual optical image. The 
residual image was obtained by subtracting the mean 2-D surface 
brightness from the original optical image (see text for details). A 
bright stellar disk is found in the nuclear region, but has no appar- 
ent association with the elongated X-ray feature. In both panels. 
North is up. East to the left. 



count s~^ arcsec"^ = (6.7 ± 0.5) erg cm~^ s~^ arcsec"^ 
(this difference is due to the decrease in sensitivity of 
the detector) ; « 80% of the unresolved counts are from 
hot gas. The total emitted luminosity from the hot gas 
is « 10^8 erg s-i in NGC 3377 and « lO^^ erg s'^ in 
NGC 4697. 

We can use these values to estimate the conversion 
factors for the other galaxies in our sample, for which 
we cannot do any spectral analysis. We can safely as- 
sume that the temperature of the hot ISM is in the 
range ~ 0.3-0.7 keV. The relative contribution of hot 
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r [arcsec] 

Fig. 8. — Count-rate surface density profiles of the unresolved 
emission in the 0.3-1.5 keV band (units of ACIS-S counts s" ^ 
arcsec"^), plotted as a function of angular distance on the sky. 
From the brightest to the faintest, symbols (and colors for the on- 
line version of the paper) are as follows: filled circles for NGC 4697 
(green for the 2000 observations, magenta for the 2003-2004 ob- 
servations); open squares (black) for NGC 3377; open circles (blue) 
for NGC 5845; open triangles (red) for NGC 4564; crosses (cyan) 
for NGC4486B. 



gas and LMXB emission varies from galaxy to galaxy 
(see a discussion in Kim & Fabbiano 2003); our choice 
of a soft energy band (0.3-1.5 keV) should reduce the 
uncertainty caused by LMXB contamination. The con- 
version factor between unresolved count rate and flux of 
the truly diffuse emission depends strongly on the detec- 
tion threshold of the point sources: we used NGC 4697 
to estimate this effect. When the point-source detection 
limit is 1 X 10'^'' erg s-^, we find that « 80% of the un- 
resolved counts in the 0.3-1.5 keV band are from the 
ISM, as mentioned earlier; « 50% for a detection limit of 
,38 o-i. g^j^j^ Qj^jy 35% for a detection limit of 



1 X 10"" erg s 

3 X 10'^^ erg s~^. In conclusion, taking into account these 
arguments, we adopt a conversion factor between count- 
rate surface density and bolometric surface brightness of 
1 count s"! arcsec"^ in the 0.3-1.5 keV band — (3.0±1.5) 
erg cm"^ arcsec"^ 
and NGC 5845. 



from the hot gas, for NGC 4564 



4.2. Density of the hot ISM 

The gas surface brightness S{r) is the projection on 
the sky of the plasma emissivity e « ngA(Tgas) where rie 
is the gas electron density and A(Tgas) the cooling func- 
tion at a gas temperature Tgag. We can obtain the gas 
density Ue by de-projecting the surface brightness, if we 
know the radial dependence of Tgas- The problem can be 
simplified by assuming that the temperature is constant, 
and that the radial profile of the surface brightness can 
be parameterized with isothermal /3-models (Cavaliere & 
Fusco-Femiano 1976, 1978) of the form: 



S{r)=So l + (r/r,) 



0.5-3/3 



(1) 



where the fitting parameters are Sq (central surface 
brightness), Tc (core radius) and /3. Physically, it can 
be shown that (3 ~ (cr^/Tgas) where a is the stellar veloc- 
ity dispersion (e.g., Ettori 2000 and references therein). 




r [pc] 

Fig. 9. — Count-rate surface density profiles of the unresolved 
emission in the 0.3-1.5 keV band (units of ACIS-S counts s" ^ 
arcsec"^), plotted as a function of projected physical distance from 
the nucleus. Symbols and colors are as in Figure 8. 



TABLE 4 

Best-fit parameters for the spectrum of the unresolved X-ray 
EMISSION IN NGC 4697. Spectral model: wabsGai*"abs*(vapec+po). 
Values in brackets were kept fixed. Errors are 90% confidence 

LEVELS for one INTERESTING PARAMETER (Ax^ = 2.7). 



Parameter 



Model 1 
Value 



Model 2 
Value 





(2.1) 


(2.1) 




(0.0) 


fj 9 + 6.6 

^■^-2.2 


rb 


(1.7) 


(1.7) 


Afpl at 1 keV 


SAtl^ X 10"'' 


8.1I22 X 10-'' 


feTgas (keV) 




r, 09 + 0.03 

""-"^-O.OS 


•^C,n(-^0) 


(1.0) 


9 1+^-3 




(1.0) 


(1.0) 


•^Pc,Ni(^0) 


(1.0) 


1 6+1° 

J^-°-0.6 


K d 

-1 »- vapcc 


3.2+°'l X 10-5 


2.7lJ;g X 10-5 


xVdof 


43.0/51 


30.7/48 


ix, 0.3-10'^ 


2.010-J X 10^9 


2.3t°| X 1039 


i'x.O.S-l.s' 


l.Sto l X 1039 


l.eto.l X 1039 




1.0l°;i X 10^9 


l.3l°;i X 1039 



^From Dickey & Lockman (1990). 
^Fixcd to the typical LMXB value of 1.7. 
""Units of ph. keV"^ cm"^ 

^Apec/vapec model normalization K^apcc — 10~"'"'^/{47r [d_4 (1 + 
•^)]^} / nf^nndV , where dA is the angular size distance to the source 
(cm), Tie is the electron density (cm~^), and uh is the hydrogen density 
(cm-^). 

■^Total unabsorbed luminosity of the unresolved emission in the 0.3—10 
keV band; units of erg s~^. 

^Total unabsorbed luminosity of the unresolved emission in the 0.3— 
1.5 keV band; units of erg s~^. 

^Unabsorbed luminosity of the diffuse thermal-plasma component in 
the 0.3-1.5 keV band; units of erg s~^. 
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Fig. 10. — Left panel: adaptively-smoothed ACIS-S image of NGC4697, in tiie 0.3-8 keV band. The ellipse marks the D25 isophote; the 
box (dashed white line) is shown in greater details on the right. Right panel: true-color X-ray image of the inner region of NGC 4697. Red 
= 0.3—1 keV; green = 1-2 keV; blue = 2—8 keV. In both panels, North is up, East to the left. 
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channel enengy (keV) 

Fig. 11. — Spectrum of the unresolved emission from NGC 4697, 
with best-fit model and residuals. The data have been fit- 
ted with a variable-abundance thermal plasma model (vapec in 
XSPEC) plus a power-law. The fit parameters are listed in Table 
4. 
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Fig. 12. — Spectrum of the unresolved emission from NGC 3377, 
with best-fit model and residuals. The data have been fitted 
with a thermal plasma model (apec in XSPEC, with fixed solar 
abundance) plus a power-law. The fit parameters are listed in 
Table 5. 



TABLE 5 

Best-fit parameters for the spectrum of the unresolved X-ray 
EMISSION IN NGC 3377. Spectral model: wabs*(apec+po). Values in 
brackets were kept fixed. Errors are 90% confidence levels for 

one interesting parameter (Ax^ = 2.7). Parameters and units 
defined as in Table 4. 



Parameter 



Value 



F 

TVpi at 1 keV 
fcTgas (keV) 
Z(Zq) 

xVdof 



(2.9) 
(1.7) 

S.ltl'l X 10-'' 

+0.19 
U.O(_Q 27 

(1.0) 

2.51^-5 ^ 

1.4/5 



ix, 0.3-10 


4.2: 


h3.2 
-2.2 


X 10^8 


ix, 0.3-1. 5 


2.0: 


hl.O 
-0.8 


X 10^* 




1.0: 


hl.4 
-0.5 


X 10^8 



^fixed to the Galactic liiic-of-sight value, units of 10 cm ; from 
Dickey & Lockman (1990). 



In this approximation, one can derive an analytical ex- 
pression for the density profile (Ettori 2000): 

2I -3/3/2 

ne(r) = no 1 + (r/rc) 



where 



-0.5 



r(3/3) 



r(3/3-0.5) reA(Tgas)' 



(2) 
(3) 



In equation 3, r(a;) is the gamma function, and the cen- 
tral surface emissivity Eq (luminosity per unit surface 
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TABLE 6 

Fitting parameters of the isothermal /3-model for the extended soft X-ray emission, and central gas densities. Errors 

ARE 90% confidence LEVELS FOR ONE INTERESTING PARAMETER (Ax^ = 2.7) 



Galaxy 


Det. limit 


Counts 

I L/ . -L . iJ JVC V 1 


I CL O aji.t_.OCv_. 




So 

(erg cm ^ s ^ arcsec ^) 


I CLi CoCC J 


(nr) 


/3 


no _ 

^ cm 3 \ 


N821'' 


2 X 10^* 








Too faint for spatial analysis 








03+0 0' 


N3377 


4 X 10^'' 


152 ± 18 


(15.9 ±3.1)10 


-6 


(3.8 ±0.8)10-1'' 


4.0 ±0.5 


215 ±25 




0.022 ± 0.004 


N4486B 


1 X 10^* 




No sig 


nificant detection of unresolved 


emission 






< 0.01 


N4564 


2 X 10^* 


25 ±8 


(7.7 ±4.2)10" 


6 


(2.1 ± 1.3)10-1' 


8.7^2-3 


630+1'" 




0.011 ±0.006 


N4697'' 


(1 X lO^'^)'^ 


1038 ± 62 


(19.2 ± 1.7)10 


-6 


(9.3 ± 1.2)10-1' 


6.1 ± 1.7 


344 ± 100 


40+°-°l 


0.017 ±0.003 


N4697<^ 


(1 X lO^'^)': 


2167 ± 92 


(14.0 ± 1.0)10 


_6 


(9.4 ± 1.2)10-1' 


4.6 ±1.7 


260 ± 100 


07+0.01 


0.019 ±0.003 


N5845 


3 X 10^** 


55± 12 


(10.7 ±3.6)10 


-C 


(3.3 ± 1.0)10-1' 


5.5±1.1 


690 ± 150 


59+010 


0.011 ±0.004 



''See Fabbiano et al. 2004. 
^2000 Jan observation. 

"^Detection limit for the point sources from all datascts combined. 
•^Combined 2003-2004 observations. 



area of the emitter, units of erg s-^ cm"^) is simply re- 
lated to the central surface brightness 5*0 (received flux 
from a unit angular element on the sky, units of erg 
s-^ cm-^ arcsec"^) by a constant geometrical factor^°, 
Eo = 5.34 X lO^^o. 

In the case of hot gas in elliptical galaxies, the gas tem- 
perature varies at most by a factor w 2 over the extent 
of the diffuse emission (e.g., Sivakoff, Sarazin & Carlin 
2004; Sivakoff et al. 2003), while the gas density can vary 
by a few orders of magnitude. Hence, taking into account 
the limited number of detected counts, we can accept the 
isothermal model as a good approximation for our case. 
We have now all the elements to determine the radial 
density profile n{r) and central density no from equa- 
tions 1-3, in the isothermal /3-model approximation. For 
each galaxy, we use the surface brightness profile inferred 
from the count rates of the unresolved emission, as de- 
scribed earlier, together with the cooling function A(r) 
from Sutherland & Dopita (1993). The results are given 
in Table 6, and the inferred radial profiles of the electron 
density (equation 2) are plotted in Figure 13. 

The above method for the deprojection of the sur- 
face brightness is based on the assumption of a spe- 
cific analytic form for the radial profile (/3-model). It 
is also possible to obtain an estimate of the density pro- 
file without assuming any models, using an onion-skin 
approach^ ^ together with the best-fit emission measure 

10 To obtain this conversion factor, recall that the area of an 
element of sky seen under a solid angle of 1 arcsec^ corresponds to 
a physical size of (48.5pc)2 x (d/10Mpc)2 = 2.35 X 10-ii(d/cm)2 
cm'^, where d is the distance to the source. The flux S received from 
a unit angular element of sky is l/(47rd^) times the luminosity emit- 
ted from the same unit angular element, or (2.35 X 10— ilci'^)/(47rd-^) 
times the luminosity E emitted from a unit surface area of the 
source. Inverting this relation gives the numerical factor between 
surface emissivity and surface brightness used in the text. 

11 An onion-skin approach consists in dividing the observed re- 
gion into concentric shells, projected as concentric annuli on the 
sky. The problem is solved from outside in. From the measured 
average temperature and luminosity of the gas in the outermost an- 
nulus, one calculates the average gas density in the outermost shell. 
The observed emission from the second annulus is given by emis- 
sion from the second shell, plus some emission from the fraction of 
the outermost shell projected in front of it; the latter component 
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Fig. 13. — Density profiles of the X-ray emitting gas, as a func- 
tion of projected physical distance from the central BH. Dotted 
(black) line: NGC 3377; dash-dotted (red) line: NGC 4564; dashed 
(green) line: NGC 4697; solid (blue) line: NGC 5845. The four 
curves have been obtained by deconvolving the best-fit /3-model 
profiles. In addition, for three of the four galaxies, we overplot 
the density directly estimated from the X-ray emission in three 
concentric shells, using an onion-skin approach, with no model as- 
sumptions. /3-models and empirical onion-skin approach give a 
consistent estimate of the hot-gas density at distances > 150 pc; 
inside that radius, the degree of flattening in a core is somewhat 
model-dependent . 



from the fitted spectral model in XSPEC (recalling that 
-f^apec ~ / nIdV). We did so for three galaxies in our 
sample^^: we defined three concentric shells, including 
the whole surface where diffuse emission is detected (i.e., 
to a radial distance of w 60"), and calculated the aver- 

is by now known, and can be subtracted. And so on, until the 
innermost annulus, whose emission comes from the innermost shell 
plus all other shells projected in front of it. 

1^ in NGC 4697, the diffuse emission extends to much larger dis- 
tances; hence this onion-skin method would give meaningful results 
only by using a larger number of shells. A more detailed analysis 
of the radial profile of the diffuse emission in NGC 4697 is beyond 
the scope of this work, which focuses on the central gas density; 
see also Sarazin et al. (2001). 
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age count rate and luminosity from each one. From the 
best-fit emission measure, we then inferred the average 
gas density in the three shells. The results are in agree- 
ment with the /3- model estimates (Figure 13), implying 
gas densities ~ 1-2 x 10~^ cm~^ at ^ 100 pc from the 
SMBH. 

Estimating the hot-gas density in the two remaining 
galaxies (NGC 821 and NGC 4486B) is more difficult, be- 
cause of the limited number of photons. For NGC 821, 



an average density n,, =4.1 



+ 10.9 
-1.5 



X 10 cm was de- 



termined in Fabbiano et al. (2004). However, this value 
was estimated over a region of radius 20", excluding the 
point sources. Therefore, it underestimates the central 
density. By scaling this value in analogy with the density 
profiles inferred for the other four galaxies, we estimate 
that it corresponds to central densities Ue = 0.03lo'g2 
cm-3 (i.e., log Tie = -1.5 ± 0.5). For NGC4486B, a cen- 
tral density w 10~^ cm~^ is the upper limit. 

For the typical values of ISM temperature and density 
in our target galaxies, the characteristic radiative cool- 
ing timcscale inside race is ~ 10^ yr, in the absence of 
external heating. In fact, both Type-la SNe and feed- 
back from the SMBH are possible heating sources, as 
we shall discuss in Section 6.2. The radiative cooling 
timescale is much longer that the accretion timcscale (^ 
free-fall timescale), which is < 10^ yr within the sphere 
of influence of the SMBH (defined as in Section 4.3). 

4.3. Bondi inflow rate of hot gas 

From the hot-gas density (determined in Section 4.2), 
we can estimate the accretion rate onto a point-like 
mass within a uniform gas, in the standard model of 
steady, spherically symmetric accretion (Bondi 1952; 
Prank, King & Raine 2002; Ruffert & Arnett 1994): 

(5-37)7(27-2) 



Mb = ttG^M^h 



cf (do) \^5 — 37 



(4) 



where 7 is the polytropic index (1 < 7 < 5/3), Cg = 
(7fcT//xmp)^/^ is the sound speed of the gas, m,, is the 
proton mass, /inip is the mean mass per particle of gas. 
We take /i = 0.62 (solar composition) and 7 = 1.33 for 
consistency with Pellegrini (2005). For the gas temper- 
ature and density "at infinity" we take their values at 
the Bondi accretion radius race = 2GMbh/Cs(oo), which 
defines the region inside which the gravitational pull of 
the SMBH dominates over the internal thermal energy of 
the gas. Henceforth, we shall refer to this region simply 
as the SMBH "sphere of influence" 

A spherically-symmetric model is probably not appli- 
cable to the behaviour of the accretion flow very close to 
the SMBH (e.g., Frank et al. 2002; Krumholz, McKce & 
Klein 2005, 2006); however, we can still take the Bondi 
rate from equation (4) as a first-order approximation of 
the rate at which the hot ISM is captured by the sphere 
of influence of the BH, at r = race- As we shall dis- 
cuss later, not all of the inflowing, hot gas reaches the 
BH, and, conversely, there are other sources of gas inside 



It is also possible to define a dyiiaiiiical sphere of influence, 
with a radius ~ GM^n/a^, where a is the stellar velocity dis- 
persion. In practice, the two definitions are very similar for this 
class of systems, because the stellar velocity dispersion ^ thermal 
velocity of virialized hot gas. 



We can more usefully express the accretion radius and 
Bondi rate as: 



Ml 



BH 



V lO^Mo 



0.5 keV 



kT 



pc. 



(5) 



Mb = 1.56 X 10" 



no 



/ M] 



BH 



0.5 keV 



VlO8M0^ 
T3) MQyr-\ 



kT 



3/2 



(6) 



VO.Ol cm- 

The numerical values for our target galaxies arc listed 
in Table 7, Cols. 3,4,5 For the hot gas phase, we used 
a temperature of 0.46 kcV for NGC 821 (Fabbiano et al. 
2004), 0.57 keV for NGC 3377, 0.33 keV for NGC 4697, 
and 0.5 keV for the other two galaxies with a detected 
hot-gas component. 

4.4. Is the central gas density underestimated? 

The hot-gas densities inferred from the deprojection 
of the diffuse X-ray emission in our target galaxies are 
comparable (though, on the lower end) with those found 
in Pellegrini's (2005) sample (ng ^ 0.01 0.5 cm^'^) and 
with the hot-gas density recently estimated for NGC 3379 
(rie ~ 0.01 cm~^: David et al. 2005). There are only two 
nearby galaxies for which the accretion radius has been 
clearly resolved from Chandra images. For the Milky 
Way (race = 1'.'8), the hot gas density ne(racc) ^ 10^ 
cm^'^ (Baganoff et al. 2003; Quataert 2003; Cuadra et 
al. 2005). For M87 (race ~ 1"9), the central density is 
« 0.17 cm~^ (Di Matteo et al. 2003). Four more galax- 
ies have accretion radii race ~ 1": in M31, a central 
density Ueir < 200 pc) ^ 0.1 cm^"^ has been estimated 
(Garcia et al. 2005; Dosaj ct al. 2002); rig « 0.04 cm'^ 
in the center of NGC 5128 (Evans et al. 2004; Kraft et 
al. 2003); rig ~ 0.1 cm~^ for the Sombrero galaxy (Pel- 
legrini ct al. 2003a); and - 0.5 cm"^ for NGC 4649 
(Soldatenkov, Vikhlinin & Pavhnsky 2003). 

Our target galaxies are slightly more distant and/or 
less massive, and, as a result, their accretion radii are 
not resolved. The observed radial profiles provide a good 
estimate of the gas distribution at distances > 150 pc; 
however, we do not have any direct measurements at 
smaller radii (Figures 9, 13), because of a lack of spatial 
resolution. If we want to determine the amount of gas 
available for accretion inside the sphere of influence of 
the SMBH, wc need to estimate the density at radii ^ 10 
pc. If /3 models are a good approximation, the presence 
of a core implies that the density stays approximately 
constant inside the inner r ~ 100 pc, making this task 
easier. But we cannot a priori rule out the possibility 
that the gas density distribution has a cusp in the nuclear 
region, for example in correspondence of a massive star 
cluster, or because cooling gas has accumulated there. 
Although we cannot resolve the nuclear gas emission in 
the Chandra images, we know that it does not domi- 
nate the point-like nuclear sources, which have harder, 
AGN-likc spectra (see, e.g., the X-ray colors in Figure 
3). A conservative constraint is that the thermal emis- 
sion inside ~ 10 pc has an X-ray luminosity < lO'^* erg 
s~^. For typical emissivities corresponding to temper- 
atures ~ 0.3-0.6 keV, this gives an order-of-magnitude 
constraint on the "true" central density rie < 10 cm~^ 
for the five target galaxies in which hot gas is found. 
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Moreover, the hot. X-ray emitting ISM may represent 
only a small fraction of the gas in the inner region. That 
can be the case in systems where gas can cool efficiently 
(cooling timescale < accretion timescale), or, vice versa, 
if gas is injected into the inner regions in a cool or warm 
phase (for example, through stellar winds) and is ac- 
creted before it has time to virialize. Cool and warm gas 
would of course elude an X-ray investigation. Therefore, 
the Bondi inflow rate of the hot-ISM component into 
the sphere of influence of the SMBH represents a firm 
lower limit to the total gas mass available for accretion. 
We shall compare the SMBH luminosity expected for gas 
accretion at the Bondi rate, in the case of a radiatively- 
efficient inflow (standard disk, Shakura & Sunyaev 1973) 
and for radiatively inefficient solutions. Conversely, we 
shall also estimate the accretion rate (expressed as a frac- 
tion of the Bondi rate) required to explain the observed 
luminosities, if we assume a radiatively-inefficient sce- 
nario. In Paper II, we shall provide a complementary 
estimate of the gas injection rate into the SMBH sphere 
of influence, based on optical studies of the stellar popu- 
lation inside the accretion radius, following the example 
of Fabbiano et al. (2004). 

5. DISCUSSION 
5.1. Accretion rates and efficiencies 

Let us define a total accretion power Pace = rj'Mc^, 
where r/' is the total accretion efficiency (radiative plus 
mechanical), and M is the rate at which matter actually 
accretes onto the BH. One immediate objective of our 
analysis is to determine whether Mb derived from X-ray 
observations of hot gas is a good tracer of M. The ef- 
ficiency rj' < 0.3, with its exact value depending on the 
physics of the accretion flow as well as on the BH spin 
parameter; more generally, we shall talk of "efficient" 
accretion when 77' ~ 0.1. We then define a bolometric 
luminosity Lboi = r/Mc'^, where rj is the radiative effi- 
ciency. In general, it is 77 = f^rj' < rj' where /r is the 
fraction of the accretion power released as radiation. In 
standard accretion models, e.g., for an accretion disk ex- 
tended down to the innermost stable orbit around the 
BH, rj w 7]' ^ 0.1; however, it is possible that 77 <C if 
most of the accretion power is released, for example, in 
a jet. Finally, only a fraction /x < 1 of the bolometric 
luminosity is emitted in the Chandra X-ray band, 0.3-10 
keV; for accreting SMBHs in typical AGN, fx ~ 10% 
(Elvis et al. 1994; Ho 1999). We shaU also define a di- 
mensionless accretion rate m = O.lMc^/LEdd) so that 

Lhoi/LEdd = (ry/0. 1)771. (7) 

Similarly, we define rfiB = 0.1AfB'^^/-^Edd- 

From the Chandra observations, the X-ray luminosity 
of the nuclear sources in our target galaxies is extremely 
sub-Eddington: Lx/^Edd = /xiboi/iEdd ^ 10~*-10"'^ 
(Table 7). From the estimated densities of the hot ISM, 
we infer dimensionless Bondi accretion rates rns ^ 10~^ 
(Table 7). 

Followig the example of previous studies (e.g., Pelle- 
grini 2005), we can compare the observed X-ray lumi- 
nosities with the inferred Bondi inflow rates. For the 
nuclei of our target galaxies, the X-ray observations im- 



ply the following constraint: 

— ^ = /xA f ^ V— 1 - 5 X 10-4 - 10- (8) 
0.1Mbc2 Vo.iy KrhBj ^ ' 

(Table 7, Col. 8). Theoretical arguments and observa- 
tional evidence in a variety of accreting galactic nuclei 
suggests that /x varies at most within a factor of w 2- 
3. The crucial problem is to constrain the other three 
parameters (/r, rj' and tti/ttib) in equation (8), to make 
them consistent with the observed SMBH luminosities. 
This will be a test for the predictions of different ac- 
cretion scenarios. We shall take into account both our 
new sample of six galaxies, and another eighteen galax- 
ies for which Bondi rates and SMBH luminosities have 
been measured or constrained (Pellegrini 2005; Garcia et 
al. 2005; David et al. 2005). 

5.2. Advective solutions (ADAF) 

Standard disk accretion predicts /r ~ 1, ~ 0.1 and 
rh/niB ~ 1. This is ruled out by the observations: if 
the gas available were efficiently accreted at the Bondi 
rate, we would expect X-ray luminosities ^ 10 100 times 
higher than observed in our target galax;ies. Moreover, 
there is a theoretical reason why we expect standard ac- 
cretion to be ruled out: in the limit of low accretion rates 
(m < 0.01), the accretion flow cannot cool efficiently 
within the infall timescale; therefore, most of the grav- 
itational energy is carried by the ions (protons and nu- 
clei) and advected into the BH, rather than being trans- 
ferred to the electrons and radiated away. This regime is 
known as Advection-Dominated Accretion Flow (ADAF, 
Narayan & Yi 1994; Ichimaru 1977). 

In the simplest, self-similar ADAF model, the luminos- 
ity Ladaf « 0.1Mc'^{rh/a'^) « {m/a)'^LEdd (Narayan & 
Yi 1994), where a is the dimensionless viscosity param- 
eter, with 0.1 < a < 1; a value a « 0.1 is typically as- 
sumed for galactic SMBHs. This corresponds to a radia- 
tive efficiency 7/ ~ OAifi/a^ ^ lOrli < 0.1. A more rigor- 
ous calculation (Merloni et al. 2003), taking into account 
synchrotron and inverse-Compton contributions, indi- 
cates that the X-ray luminosity ^x.adaf oc w^'^L^I^. 
In the parameter range of our study, the two estimates 
agree within a factor of 3. Moreover, in the ADAF sce- 
nario, the radial inflow velocity inside the Bondi radius is 
only ~ acg. For a fixed hot-gas density and temperature 
at infinity, this implies that a more accurate estimate of 
the accretion rate is 7n ~ a7hB ~ 0.l7nB (Narayan 2002). 

Using the ADAF values of efficiency and accretion rate 
(from suitable rescalings of the Merloni et al. 2003 model 
fit), we estimated the 0.3-10 keV luminosity predicted 
by the advective model for our target galaxies (Table 7, 
Col. 9). Conversely, we also determined the accretion 
rates that would be required to reproduce the observed 
luminosities, at an assumed ADAF-like efficiency (Table 
7, Col. 10). We find (Table 7 and Figure 14) that, at an 
accretion rate 777 « 0.1777b (as predicted by the ADAF 
model), the SMBH hmiinositics arc underestimated by 
a factor of 10 (for NGC4697) and by three to four or- 
ders of magnitude for the other five galaxies. Conversely, 
we find that accretion rates m w 0.6-6rhB are required 
to reproduce the observed luminosities, at an ADAF-like 
efficiency. In particular, an accretion rate « O.Sttt-b suf- 
fices to explain the X-ray luminosity of NGC 4697, but 
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TABLE 7 

Physical properties of the nuclear BHs (Lx is the unaborbed luminosity in the 0.3-10 keV band) 



Galaxy 
(1) 



A^BH 
(IO8M0) 
(2) 



(arcsec/pc) (Ms yr -•-) 



(3) 



(4) 



(5) 



logLx 
(erg s-i) 
(6) 



(7) 



Lx/(0.1A^Bc2) 

~ ■f^ ((Tl) ( riiB ) 

(8) 



logi^X.ADAF M/Mb 

(M = 0.1Mb) 

(9) (10) 



N821 

N3377 

N4486B 

N4564 
N4697 
N5845 



0.85 
1.0 



+0.36 
0.35 
+0.9 
-0.1 



[^6.0 
0.5 

0.56 



+3.0 
2.0 
+0.5 
0.2 



+0.03 
0.08 



1 7+0.2 
^•'-0.1 

2 4+°-" 



0'.'07/7.8 
0^14/7.4 

0^'61/50.4 
0('05/4.2 
0^'06/4.7 
0'.'38/21.6 
OC 16/20.2 



3.8 X 10-5 

2.8 X 10-5 

< 1 X 10-^ 

< 1 X 10-5 

5.4 X 10-6 

1.5 X 10-4 

9.9 X 10-5 



-4 70+105 



-4.92 



< 



+0.50 
-0.24 

4.1 



< -5.0 
-5.38 



+0.40 
0.40 



-4 41+0.11 
^•^■■^-o.io 

-4 74+0.37 
^•'*-0.44 



< 38.7 
38.5 
38.4 
38.4 
38.9 
38.6 
39.4 



< -7.3 
- -7+0.10 

'•^'-0.39 



-8.5 
-7.44; 
-6.96 



+0.2 
0.3 

+0.28 
0.39 
+0.23 
0.22 



7 74+0.03 
'•'*-0.06 

-7 00+0-26 



< 6 X 10-3 
2.3 X 10-3 

> 4 X 10-5 

> 4 X 10-3 

0.03 
4.6 X 10-4 
5.6 X 10-3 



36.4 
36.0 

< 38.2 

< 35.2 
34.3 
37.5 
36.8 



< 2 
1.3 

> 0.1] 
> 2 
6 

0.3 
1.4 



Note. — CoL(l): galaxy ID. NGC4486B is listed twice: the first line (in square brackets) assumes the mass of Kormendy et al. (1997); the 
second line assumes the lower mass indirectly derived in our Paper II. Col. (2); SMBH mass, see Table 1. Col. (3): accretion radius, from equation 5. 
Col. (4), (5): Bondi inflow rates estimated from the density and temperature of the diffuse hot ISM. Col. (6), (7): unabsorbed X-ray luminosity of 
the nuclear source, in the 0.3—10 keV band. Col. (8): comparison between observationally-dctcrmined X-ray luminosity and bolometric luminosity 
expected from standard, radiatively-cfflcicnt accretion. Col. (9): prediction for the 0.3—10 keV luminosity for an ADAF accretion flow, with a 
viscosity parameter a — 0.1; values were suitably rescaled from the flt of Merloni et al. (2003), using an accretion rate rh Si ams, as expected for 
an ADAF scenario. Col. (10): accretion rates required to satisfy the observed (unabsorbed) X-ray luminosities, as a fraction of the Bondi inflow 
rate of hot gas. 



rates m > rfiQ arc required for the other five galaxies. 
If tlie liot-ISM inflow was the only source of gas, such a 
high accretion rate would be difficult to reconcile with 
the ADAF prediction. We take this finding as a strong 
indication that, for our target galaxies, the gas avail- 
able for accretion has been underestimated: there are 
other sources of fuel for the SMBH, such as stellar winds, 
that could be more significant than the Bondi infiow; this 
would explain why we infer to > tob- We shall discuss 
this scenario in detail in Paper II (see also Fabbiano et 
al. 2004; Pellegrini 2005). 

More generally, there are « 20 quiescent or very faint 
SMBHs for which the SMBH mass and X-ray luminos- 
ity, and the density of the surrounding hot ISM, have 
been estimated or constrained (Pellegrini 2005; Garcia 
et al. 2005; David et al. 2005). When plotted in an X- 
ray-luminosity versus Bondi-rate diagram (Figure 14), 
the distribution of the datapoints does not provide any 
obvious correlations between these two parameters, or 
between the luminosities observed and predicted by the 
ADAF model. (The standard disk model fares much 
worse, and is clearly ruled out for all but one galaxy). 
Some SMBHs, including those in the six galaxies of our 
sample, tend to be more X-ray luminous than predicted 
by the standard ADAF model (i.e., they lie above the 
predicted to « O.Itob ADAF line in Figure 14). This 
seems to be the case for most galaxies with extremely 
low hot-ISM densities in the nuclear region (rig ^0.1 
cm-^), and therefore very low Bondi accretion rates 
(rhe < 10~*). On the other hand, SMBHs in galaxies 
with estimated Bondi accretion rates tob ^ IQ-^ tend to 
lie below the ADAF line (Figure 14); that is, the ADAF 
model overpredicts their observed luminosities, in some 
cases by more than two orders of magnitude. 

In the framework of inefficient accretion, we can ex- 
plain the former group of galaxies by considering addi- 
tional sources of gas which boost the accretion rate above 



the classical Bondi value, as we shall discuss in Paper II. 
Conversely, an additional mechanism must be at work in 
the latter group of galaxies, reducing the true accretion 
rate to values rh < IQ-^tob- As an aside, we note that 
the lack of SMBHs in the bottom left of the diagram 
(i.e., with m < 10^'''tob and Lx < lO^^^iEdd) rnay be 
partly due to a selection effect, because both the nuclear 
point-like source and the diffuse hot-ISM emission would 
be close to or below the Chandra detection limit in those 
galaxies. 

5.3. Convection and outflows 

The original ADAF solution has subsequently been de- 
veloped by taking into account two physical effects that 
are likely to further reduce the accretion rate, so that 
m ^ Tob (see Narayan 2002 and references therein, for 
a review). Firstly, the infiowing gas is likely to become 
convectively unstable due to its inefficient cooling (Con- 
vection Dominated Accretion Flow, CDAF: Narayan et 
al. 2000); convective instabilities are also likely to oc- 
cur when a magnetic field is frozen into the infalling gas 
(Igumenshchev & Narayan 2002); or when the inflow- 
ing gas has a large angular momentum (Proga & Begel- 
man 2003) or vorticity (Krumholz, McKee & Klein 2005, 
2006). Secondly, inefficient cooling causes the gas to 
have positive Bernoulli parameter, thus triggering large- 
scale winds and outflows which reduce the amount of gas 
reaching the BH (Advection Dominated Inflow-Outflow 
Solutions, ADIOS: Blandford & Begelman 1999). 

Low, ADAF-like efficiencies and a reduced accretion 
rate with respect to the Bondi rate (to/tob ^ 10-^- 
10^^), via CDAF and/or ADIOS solutions, is a possible 
explanation for those galaxies in which hot-gas densities 
are ~ 0.1-10 cm-^, and correspondingly, IQ-^ < me ^ 
10~^. They lie below the standard ADAF line in Figure 
14. Among them, the SMBH in the Milky Way, for which 
Tob « 2 X 10^^; an accretion rate to/tob ^ 10^"^ has been 
suggested (Quataert 2003). The SMBH in the nucleus 
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Fig. 14. — Top panel: relation between normalized X-ray lu- 
minosity Lx/^Edd E'lid normalized Bondi accretion rate rhe = 
A/e/A/Edd- Bottom panel: zoomed-in view of the same diagram, 
to allow a better identification of the six galaxies of the present 
study (datapoints plotted in red in the online version). The data- 
points in black have been collected by Pellegrini (2005) (see refer- 
ences therein for the original measurements; a slightly revised mass 
has been adopted for NGC4636, from Wang & Zhang 2003). The 
large filled circles (in blue, in the online version) represent Sag A* 
(Baganoff et al. 2003) and M31* (Garcia et al. 2005). The data- 
point plotted in green is NGC 3379 (David et al. 2005). We have 
compared the observed X-ray luminosities with those predicted, as 
a function of Bondi accretion rate, by: a) standard radiatively- 
efficient accretion {tj = 0.1, rh = ifi^: solid line); b) standard 
ADAF solution (rj = lOrh, rh = O.lrhe: dashed-dotted line); c) 
a scenario with ADAF-like efficiency {rj = 10m) but higher accre- 
tion rate (m = ttib: dashed line). SMBHs with X-ray luminosities 
falling below the standard ADAF line are likely to require an ad- 
ditional mechanism to reduce the true mass accretion rate further 
below the Bondi value (e.g., ADIOS or CDAF solutions). Con- 
versely, for the six target galaxies studied in this paper, the Bondi 
rate is too low to be consistent with the observed X-ray luminosi- 
ties, if the accretion process has ADAF-like radiative efficiencies. 



of M31 also belongs to this group (Garcia et al. 2005). 
At higher central gas densities (> 10 cni~^, correspond- 
ing to tob ^ 10~^), theoretical and observational argu- 
ments show that accretion will become radiatively effi- 
cient, with the formation of a standard, optically-thick, 
geometrically-thin accretion disk. SMBHs in this regime 
will appear as fully- fledged AGN. 



5.4. Intermittent accretion 

We have assumed so far that the accretion rate can be 
expressed as a constant fraction of the gas available. We 
also need to take into account the possibility that accre- 
tion is intermittent, so that the system may go through 
phases when the SMBH accretion rate is much lower than 
the Bondi inflow rate, and subsequent phases in which it 
can be much higher. The nuclear region would alterna- 
tively build up or be depleted of gas. 

In one version of this scenario, it is the feedback from 
the nuclear X-ray source that drives and self-regulates 
the activity cycles: when the nuclear source is active, 
the surrounding ISM is radiatively heated (Ciotti & Os- 
triker 2001; Ostriker & Ciotti 2005), or, in the case of 
radiatively-inefficient accretion, is heated by the bulk ki- 
netic energy of jets, outflows or convective motions (Bin- 
ney & Tabor 1995; Omma et al. 2004), until the radial 
inflow is stopped and the nucleus switches off. Thus, the 
Bondi inflow estimated from the gas density at large radii 
(at the limit of the Chandra resolution) may not provide 
a good indication for the time-dependent accretion rate 
and SMBH luminosity. This is a possible explanation for 
the observed lack of any relation between ix, -^^bh and 
Mb (Pellegrini 2005). 

Alternatively, accretion may cycle between radiatively 
inefficient flows (e.g., CDAF) in which gas accumulates 
inside the accretion radius, and a standard, radiatively- 
efficient accretion disk in which the stored gas is quickly 
accreted into the BH. However, this transition should 
only happen at accretion rates m 10^^, much higher 
than estimated in our target galaxies. 

Intermittent accretion may also be due to thermal- 
viscous ionization instabilities, triggered by an ionization 
front propagating inwards or outwards across an accre- 
tion disk (Siemiginowska, Czerny & Kostyunin 1996); a 
change of the viscosity parameter by a factor of 4 between 
the two states may lead to changes in luminosity by four 
order of magnitudes. Again, this scenario requires that 
the inflowing gas can cool efficiently to form an accretion 
disk. 

It was suggested (Fabbiano et al. 2004) that intermit- 
tent ejections may be responsible for the knotty jet-like 
feature in the nuclear region of NGC 821. If this emission 
is connected with nuclear activity, the most likely expla- 
nation would be that of thermal emission from shocks in 
the hot ISM resulting from a past nuclear outburst. 

6. CONCLUSIONS 

We have used Chandra^ supported by archival HST and 
ground-based optical images, to study the nature, accre- 
tion power and mass budget of X-ray-faint SMBHs in el- 
liptical/SO galaxies. For this purpose, we have chosen a 
sample of six nearby galaxies with quiescent SMBHs, five 
of which with accurate BH mass determinations (from 
stellar kinematics). From the X-ray images, we found 
different morphologies for the X-ray emission in the nu- 
clear region. In NGC 821, a jet-like feature has been pre- 
viously suggested (Fabbiano et al. 2004). In NGC 4697, 
the nuclear emission is resolved into a few point-like 
sources (one of which is coincident with the nucleus). In 
other cases, (NGC 3377, NGC 4564, NGC 5845) the nu- 
clear emission is extended on scales of a few hundred pc 
but its morphology remains unclear. Obscuration of the 
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nuclear source may play a role in some cases; for exam- 
ple, in at least one galaxy (NGC5845), the SMBH may 
be partly obscured by an edge-on disk of dust and stars, 
and the observed extended X-ray emission may be due 
to scattering in the surrounding plasma above and below 
the disk plane. 

In five of the six galaxies, the X-ray colors of the nu- 
clear emission are hard, consistent with an AGN spec- 
trum. Hence, it is likely that most of the nuclear emis- 
sion comes from the accreting SMBH, rather than from 
a central concentration of hot gas. In addition, softer 
X-ray emission from diffuse hot gas is detected around 
the nuclear region, out to 1 kpc. We have fitted the 
soft X-ray brightness profiles with /3 models to estimate 
central hot gas densities (typical values are ^ 0.02 ±0.01 
cm~^). For two of our galaxies, the flux from the dif- 
fuse component is high enough to allow also for spectral 
fitting. NGC 4486B is the only case where the central, 
point-like X-ray source is soft but there is no evidence of 
diffuse thermal emission around it. 

The Bondi accretion rate gives an estimate of the rate 
at which diffuse hot gas flows into the sphere of influ- 
ence of the SMBH (characteristic radius race ~ 10 pc 
^ O'.'l). Typical Bondi inflow rates vary from ~ a few 
xlO~^ to ^ 1O^'*M0 yr~^ in our target galaxies This 
is a lower limit to the gas injection rate into the sphere 
of influence, because it does not take into account addi- 
tional gas sources directly inside race (for example, stellar 
winds or gas losses from Type- la supernovae). 

Based on the X-ray data, we tried to answer the fol- 
lowing (inter-related) questions: whether the observed 
X-ray SMBH luminosities are consistent with standard 
or radiatively-inefhcient accretion; whether the X-ray- 
inferred Bondi rate is a good estimate of the true accre- 
tion rate onto the SMBH; and whether there is a corre- 
lation between Bondi rate and X-ray luminosity of the 
SMBH. We found that standard accretion is ruled out for 
all galaxies, as it would overestimate the SMBH luminosi- 
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ties. On the other hand, our quiescent sample of SMBHs 
are also brighter than what is predicted by radiatively- 
inefficient solutions (ADAF, ADIOS, CDAF); we inter- 
pret this as evidence that the Bondi inflow rate of the 
hot ISM underestimates the total gas injection rate into 
the sphere of influence of the SMBH. However, when we 
consider the other X-ray faint SMBHs for which the nu- 
clear luminosity and the Bondi rate are known from the 
literature, we see that most of them are consistent with, 
or fainter than predicted by the ADAF model (as noted 
in Pellegrini 2005). In those cases, it appears that the 
Bondi inflow rate of hot gas is an over-estimate of the 
true accretion rate. Overall, there seems to be no corre- 
lation between the Bondi rate and the X-ray luminosity 
of the SMBH. A possible explanation suggested (Pelle- 
grini 2005) for this lack of correlation is that accretion is 
intermittent, self-regulated by feedback from the BH: it 
may switch between a state in which m <^ me (when gas 
builds up in the inner few pc), and one in which m 3> me 
(as it is depleted). 

In Paper II we will use optical observations to inves- 
tigate additional sources of fuel for the SMBH; this will 
provide a better estimate of the total accretion rate, and 
will allow us to re-analyze its correlation with the X- 
ray luminosity. We shall then be able to estimate the 
energy and mass budget, discuss the conditions for a 
steady-state, and distinguish between different varieties 
of radiatively-incfflcicnt models. 
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